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Time delayed feedback control is one of the most successful methods to discover 
dynamically unstable features of a dynamical system in an experiment. This 
approach feeds back only terms that depend on the difference between the current 
output and the output from a fixed time T ago. Thus, any periodic orbit of period T 
in the feedback controlled system is also a periodic orbit of the uncontrolled system, 
independent of any modelling assumptions. 

It has been an open problem whether this approach can be successful in gen¬ 
eral, that is, under genericity conditions similar to those in linear control theory 
(controllability), or if there are fundamental restrictions to time-delayed feedback 
control. We show that there are no restrictions in principle. This paper proves 
the following: for every periodic orbit satisfying a genericity condition slightly 
stronger than classical linear controllability, one can find control gains that stabilise 
this orbit with extended time-delayed feedback control. 

While the paper's techniques are based on linear stability analysis, they ex¬ 
ploit the specific properties of linearisations near autonomous periodic orbits in 
nonlinear systems, and are, thus, mostly relevant for the analysis of nonlinear 
experiments. 

1. Introduction 

Time-delayed feedback control was originally proposed by Pyragas in 1992 as a tool for 
discovery of unstable periodic orbits (one frequent building block in nonlinear systems with 
chaotic dynamics or multiple attractors) in experimental nonlinear dynamical systems [1]. 
Pyragas proposed that one take the output x(t) <G IR n of a dynamical system and feed back in 
real time the difference between this output and the output time T ago into an input u(t) G IR' lu 
of the system (multiplied by some control gams K T e R nuXn ): 

u(t) = K T [x(t —T) — x(t)]. (1) 

In a first experimental demonstration, Pyragas and Tamasevicius successfully identified and 
stabilised an unstable periodic orbit in a chaotic electrical circuit [2]. Socolar et al in 1994 [3] 
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introduced a generalisation of time-delayed feedback (which is often used in place of (1) and 
is implemented as shown in Figure 1 as a block diagram): 

u(t) = K T [x(t) -x(t)], where 
x(t) = (1 — e)x(t — T) + ex(t — T), 


and e G (0,1], called extended time-delayed feedback. If e = 1, feedback law (2) reduces to time- 
delayed feedback (1), if e = 0 feedback law (2) degenerates to classical linear feedback with a 
fixed T-periodic reference signal x(t) (see below (3) for a discussion). Note that, for example in 
[3], the variable x(t) was eliminated in the mathematical discussion by writing 

While this would suggest that knowledge of all history of x is required to initialise the system, 
in the experiment the feedback control was implemented as shown in the block diagram in 
Figure 1, which is equivalent to (2). By construction of the feedback laws (1) and (2), for 


u(t) = K 1 


Jjl-erVt-jT) 


j=i 



Figure 1: Block diagram for extended time-delayed feedback (2), as applied to an experiment, for 
example, in [3]. We prove generic stabilisablity for the case when the output x is the whole internal 
state of the dynamical system and the input u is scalar. Triangle block symbols are multiplications of 
the signal by the factor in the block. 

e > 0 every periodic orbit of period T of the dynamical system with feedback control is also 
a periodic orbit of the uncontrolled system (u = 0). 1 Flowever, the stability of the periodic 
orbit may change from unstable without control to asymptotically stable with control for 
appropriately chosen gains K. 

The delayed terms x(t — T) and x(t — T) make extended time-delayed feedback control 
different from the classical linear feedback control, which has the form 

u(t) = K T [x*(t) — x(t)], (3) 

where x* (t) is, for example, a known unstable periodic orbit of the dynamical system governing 
x. While the goal of (3) is to stabilise a known reference output (in this case a periodic orbit), 
time-delayed feedback is able to stabilise and, thus, find a-priori unknown periodic orbits. For 
this reason time-delayed feedback originated, and has found most attention, in the physics and 
science community, rather than in the control engineering community. It can be used to discover 
features of nonlinear dynamical systems inaccessible in conventional experiments, such as 
unstable equilibria, periodic orbits and their bifurcations, non-invasively. A few examples 
where time-delayed feedback (or its extended version) have been successfully used are: control 

1 For £ ^ 0, (2) with T-periodic x implies that x = x for all t. 
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of chemical turbulence [4], all-optical control of unstable steady states and self-pulsations 
in semiconductor lasers [5, 6, 7], control of neural synchrony [8, 9, 10], control of the Taylor- 
Couette flow [11], atomic force microscopy [12] and (with further modifications) systematic 
bifurcation analysis in mechanical experiments in mechanical engineering [13,14,15]. 

One difficulty for time-delayed feedback is that there are until now no general statements 
guaranteeing the existence of stabilising control gains K under some genericity condition on 
the dynamical system governing x and its input u, such as controllability. This is in contrast to 
the situation for classical linear feedback control (3), where the following is known [16]: if the 
periodic orbit x* is linearly controllable by input u in p periods (this is a genericity condition) 
then one can assign its period-pT monodromy matrix to any matrix with positive determinant 
by pT-periodic feedback gains K(t) T € R nuXn . 

The greater level of difficulty for (extended) time-delayed feedback is unsurprising since the 
feedback-controlled system acquires memory. Let us assume that the measured quantity x is 
governed by an ordinary differential equation (ODE) x.(t) = f (x(t), u(t)) (which is autonomous 
without control (u = 0) and non-autonomous with classical feedback control (3)). Then x and x 
will be governed by a delay differential equation (DDE) if u is given by time-delayed feedback 
(1), or by an ODE coupled to a difference equation if u is given by extended time-delayed 
feedback (2) with e e (0,1) (we will refer to both cases simply as DDEs). This means that the 
initial value for both, x and x, is a history segment, a function on [—T, 0] with values in R n . 
In DDEs periodic orbits have infinitely many Floquet multipliers. 2 Section 2 will review the 
development of analysis for the time-delayed feedback laws (1) and (2). This paper proves a 
first simple generic stabilisability result for extended time-delayed feedback control (2) with 
time-periodic gains K(t) (similar to results for classical linear feedback control). 

Main result The following theorem states that the classical approach to periodic feedback gain 
design by Brunovsky [17] can be applied to make (2) stable in the limit of small e > 0 in the 
simplest and most common case of a scalar input u (thus, n u = 1) and linear controllability of 
the periodic orbit by an input at a single time instant. 

Theorem 1.1 (Generic stabilisability with extended time-delayed feedback) 

Assume that the dynamical system 

x(t) = f(x(t),u(t)) (f : R n x R i-> R n smooth) (4) 

with u = 0 has a periodic orbit x* (t) of period T > 0, and assume that the monodromy 
matrix 3 Po of x* from time 0 to T is controllable with bo = 3 u f(x*(0),0) (that is, 
det[b 0 , P 0 b 0 ,..., Pj' _1 b 0 ] / 0). 

Then there exist gains Ko € R n such that x* as a periodic orbit of the feedback 
controlled system (4) with (see below for the definition of the function A§) 

u(t) = A 5 (t)Ko [x(t) -x(t)], x(t) = (1 - e)x(t-T) + ex(t-T) (5) 

has one simple Floquet multiplier at 1 and all other Floquet multipliers inside the unit 
circle for all sufficiently small e and 6. 


2 Floquet multipliers are the eigenvalues of the linearisation of the time-T map along the periodic orbit. 
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The function A§ is zero except for a short interval of length 6 every period T such that the 
feedback u has the form of a short but large near-impulse: 4 


Ag(t) 


1/5 if t|mod[0,T) £ [O' 

o if tlmod[0,T) ^ [0, 5]. 


( 6 ) 



Figure 2: Illustration of Floquet multiplier spectrum for extended time-delayed feedback with single 
input. Using the appropriate control gains Ko, n Floquet multipliers can be freely assigned up to 
determinant restrictions (n = 2 in the illustrated case). The other Floquet multipliers lie on a circle 
of radius e/2 around 1 — e/2, accumulating at 1 — e. This spectrum is achieved asymptotically for 
sufficiently short and strong impulses (6 <C 1) and small e. A simple trivial multiplier at 1 is always 
present. 


Remarks — Constant gains The gains as constructed are periodic. This is to be expected since 
there are no general results for constant gains K e lR n for the classical linear feedback case 
(3), either. Furthermore, simple examples show that the above statement can definitely not be 
made when we restrict ourselves to constant gains in (5): u(t) = K T [x(t) — x(t)] with K £ JR n . 
See Section 5 for an example. 

Properties of the spectrum of the linearisation (see also Figure 2 for an illustration) The claim of 
Theorem 1.1 is about linear stability of the periodic orbit (x(t),x(t)) = (x*(t),x*(t)) of (4)-(5). 
Thus, we have to consider the problem (4)-(5), linearised in (x(t),x(t)) = (x*(t),x*(t)): 

x(t) = A(t)x(t) + b(t)A§(t)Kg [x(t) -x(t)], 
x(t) = (1 — e)x(t — T) + ex(t — T), 

where A(t) = 3 x f(x*(t),0) and b(t) = 3 u f(x*(t),0). 

The gains Ko are identical to those chosen by Brunovsky [17] for the classical feedback 
spectrum assignment problem (note that Brunovsky made weaker assumptions on A(t) and 
b(t) than Theorem 1.1). One can choose the gains Ko to place the n Floquet multipliers 
(k = 1,... ,n) of 

x = A(t)x — b(t)A5(t)Ko x, 

4 The notation t| mo d[o,T) refers to the number t e [0, T) such that (t — t)/T is an integer. 
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anywhere inside the unit circle subject to the restriction that they have to be eigenvalues of a 
real matrix with positive determinant. 

However, DDEs such as (7) may have infinitely many Floquet multipliers. Theorem 1.1 rests 
on a perturbation argument for small e > 0 for the other, delay-induced, Floquet multipliers . 5 
At e = 0 the difference equation for x. in (7) simplifies to x(t) = x(t — T). Thus, an arbitrary 
initial history x with period T will not change under the time-T map of (7). This results for 
£ = 0 in a spectrum of (7) consisting of 

• the finitely many assigned Floquet multipliers A^ (k = 1,..., n) as determined by the 
gains Ko, and (assuming all A k / 1) 

• the spectral point A^ = 1 with an infinite-dimensional eigenspace. Specifically, if we 
choose the space of continuous functions C([—T, 0]; lR n x R n ) as phase space for (7) then, 
for £ = 0 , the eigenspace for Aoo = 1 is 

{(x,x) G C([—T,0];R n x R n ) : x(0) = x(-T), x(0) = x(-T), 

x = A(t)x + b(t)A 6 (t)Kj [x(t) -x(t)]|. 

Note that, since A^ / 1 for k = 1,.. ,,n, the ODE x = A(t)x + b(t)A 6 (t)KQ [x(tj — x(t)] 
has a unique periodic solution x for all periodic functions x. This means that for every 
T-periodic x there is an eigenvector for A<x, = 1 with this x-component. 

The general theory for DDEs [18] ensures that for positive (small) £ the Floquet multipliers 
Aic (k = 1,..., n) are only slightly perturbed, and that the infinitely many Floquet multipliers 
emerging from A^ accumulate to the spectrum of the essential part, the difference equation in 
(7) with the x terms only: x(t) = (1 — £)x(t — T). Specifically, the only accumulation point of 
the spectrum of (7) for £ G (0,1) is at 1 — £ and the stability of (7) is determined by the location 
of the Floquet mulitpliers emerging from the perturbation of A^ (of which at most finitely 
many can lie outside the unit circle). The detailed analysis in Section 3 will show that for 
small £ > 0 the Floquet multipliers emerging from A^ lie close to a circle of radius £/2 around 
1 — e/ 2, inside the unit circle (except for the unit Floquet multiplier), as shown in Figure 2 for 
n = 2. 

Trivial multiplier The eigenvector to the trivial multiplier 1 is x* (0), corresponding to the lin¬ 
earised phase shift (for every s 6 R,t i-r x*(t + s)is also a solution of the system with extended 
time-delayed feedback (4), (5)). Section 4 gives a modification of Theorem 1.1 with a function 
A 5 depending on x(t) instead of t switching the gains on and off. Then the feedback controlled 
system becomes autonomous. In this modified system with autonomous (but nonlinear) ex¬ 
tended time-delayed feedback the periodic orbit x* is asymptotically stable in the classical 
sense. 

Timing of impulse In (6) we chose the timing of the impulse (the part of the period [0, T) where 
A§ is non-zero) as [0,6] without loss of generality. The genericity condition in its most general 
form requires that there must be a time t G [0, T] such that the monodromy matrix from t to 
t + T and 3 u f (x* (t), 0) are controllable. As the uncontrolled system is autonomous, we can 
shift the phase of the periodic orbit x*, considering x*(t + •) instead of x*. 

5 The perturbation is not a small-delay perturbation since the delay T and one coefficient in front of the delay, 1 — e, 
are not small. 
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Practical considerations The result gives precise control over the Floquet multipliers in the limit 
of small 5 and e. For small 6 the feedback control corresponds to a sharp kick once per period, 
which is not practical for strongly unstable periodic orbits. Flowever, the gains found with 
the help of Theorem 1.1 provide a feasible starting point for optimisation-based spectrum 
assignment methods (continuous pole placement) as constructed by Michiels et al [19, 20] and 
adapted to time-delayed feedback (1) [21, 22, 23]. In the context of continuation one can 
combine the gains provided by Theorem 1.1 as starting points, continuous pole placement, 
and the automatic adjustment of the time delay T demonstrated in [24, 25] to create a feedback 
control that non-invasively tracks a family of periodic orbits in a system parameter. 

2. Review: analysis of (extended) time-delayed feedback 

The initial proposals of time-delayed feedback (1) and its extended version (2) were accompa¬ 
nied with demonstrations in simulations and experiments, showing that this type of feedback 
control can be successful [1, 2, 3], but not with general necessary or sufficient conditions for 
applicability or with constructive ways to design the feedback gains. 

However, it was quickly recognised that time-delayed feedback can be applied to periodic 
orbits that are weakly unstable due to a period doubling bifurcation or torus bifurcation 
[26,27], Hence, time-delayed feedback is often associated with control of chaos, because it can 
be used to suppress period doubling cascades. However, general sufficient criteria were rather 
restrictive [28], requiring full access to the state (x governed by x(t) = f(x(t)) + u(t) with 
u e ]R n ). A first general result was negative, the so-called odd number limitation for periodically 
forced systems [29], showing that extended time-delayed feedback cannot stabilise periodic 
orbits in periodically forced systems with an odd number of Floquet multipliers A with Re A > 1 
(and no Floquet multiplier at 1). This theoretical limitation is not a severe restriction in practice 
since one can extend the uncontrolled system with an artificial unstable degree of freedom 
before applying time-delayed feedback [30]. Fiedler et al showed that this limitation does not 
apply to autonomous periodic orbits [31, 32], Since then general results have been proven for 
weakly unstable periodic orbits with a Floquet multiplier close to 1 (but larger than 1, [33]), or 
near subcritical Hopf bifurcations [34, 35]. A review of developments up to 2010 is given in 
[36]. 

An extension of the odd number limitation to autonomous periodic orbits (with trivial 
Floquet multiplier) was given by Hooton & Amann [37, 38] for both, time-delayed feedback 
(1) and its extension (2). However, these limitations merely impose restrictions on the gains K. 
They do not rule out feedback stabilisability a priori (which is in contrast to the statements 
about periodic orbits in forced systems). 

3. Spectrum of linearisation for extended time-delayed feedback-controlled 
system 

Let us consider a feedback controlled dynamical system with extended time-delayed feedback 
control and arbitrary time-dependent gains K(t) 6 R n : 

*(t) = f(x(t),u(t)), (8) 

u(t) = K(t) T [x(t) — x(t)], (9) 
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x(t) = (1 — e)x(t — T) + ex(t — T). 


(10) 


This system is governed by an ordinary differential equation (ODE) without control (u = 0) 
and a delay differential equation (DDE) with control. We assume that the uncontrolled system 
x(t) = f(x(t),0) has a periodic orbit x* of period T. This periodic orbit x* is also a periodic 
orbit of (8)-(10) if e > 0: x(t) = x(t) = x*(t). System (8)-(10) is a DDE with the phase space 

{(x, x) € C([—T, 0];R n x lR n ) : x(0) = (1 - e)x(-T) + ex(-T)}. 


Hale & Verduyn-Lunel [18] treated DDEs of the type of system (8)-(10) (which contains 
difference equations) as part of their discussion of neutral DDEs. The essential part of the 
semiflow generated by ( 8 )-( 10 ) is governed by the part of ( 10 ) containing x.: x(t) = (1 — e)x.(t — 
T), which is linear and has spectral radius 1 — e. Thus, it fits into the scope of the theory as 
described in the textbook by Hale & Verduyn-Lunel [18]. Specifically, the asymptotic stability 
of the periodic orbit given by x(t) = x(t) = x* (t) is determined by the point spectrum of the 
linearisation of (8)-(10). Hence, the periodic orbit x* is stable if all Floquet multipliers of the 
linearisation along x* except the trivial multiplier 1 are inside the unit circle (and the trivial 
Floquet multiplier 1 is simple). We denote the monodromy matrix 6 of 

x = [A(t) - qb(t)K(t) T ]x(t), where A(t) = 3 x f(x*(t),0), b(t) = 3 u f(x*(t),0) (11) 

for q E C by P(q). Thus, the monodromy matrix of the uncontrolled system x(t) = A(t)x(t) 
equals P(0), which we denote by 

Po = P(0). ( 12 ) 


With this definition of P(q), Floquet multipliers of the linearisation of (8)-(10) in x* different 
from 1 — e are given as roots of 


h(A; e) := det 


AI — P 1 — 


A —(1 —e) 


(I is the identity matrix; see Section A.l for detailed proof). The following lemma states that 
the gains K(t) can only stabilise a periodic orbit x* with extended time-delayed feedback 
and small e, if they are stabilising with classical linear feedback (that is, when replacing the 
recursively determined signal xby the target orbit x*: u(t) = K(t)[x*(t) — x(t)]). (Recall that 
A(t) = 3 x f(x*(t),0), b(t) = 3 u f(x*(t),0).) 

Lemma 3.1 (Extended time-delayed feedback stabilisation implies classical stabilisation) 

If the linear system 

x(t) = [A(t)-b(t)K(t) T ]x(t) (13) 

has at least one Floquet multiplier outside the unit circle, then there exists a e max G (0,1) 
such that the periodic orbit x* is unstable for the extended time-delayed feedback 
( 8 )-( 10 ) for all e £ (0, e ma x)- 


6 Thus, P(p) is defined as the solution y at time T of the linear differential equation y (t) = [A(t) — qb(t)K(t) T ]y (t) 
with initial value y (0) = I (I is the identity matrix). 
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Proof The Floquet multipliers of (13) are given as roots of h(A;0) = det(AI — P(l)). We denote 
the root with modulus greater than 1 by Ao such that h(Ao; 0) = 0. Consequently, for all A in the 
ball B t (Ao), where r = (|Aol — 1)/2, the difference h(A; e) — h(A;0) is uniformly bounded and 
analytic for all e e (0, 1) and all A in B r (Ao). Since Ao must have finite multiplicity as a root of 
h(-; 0), h(A; e) must have a root in B r (Ao) for sufficiently small e > 0 (say, e <E (0, e max )), too. By 
choice of r this root lies outside of the unit circle. ( This ends the proof of Lemma 3.1.) □ 

Lemma 3.1 shows that gains K(t) that stabilise with extended time-delayed feedback with 
small e also have to feedback-stabilise in the classical sense. Since K is an arbitrary periodic 
function there are many ways to construct gains for the classical linear feedback control 
u(t) = K(t)[x* (t) — x(t)] for periodic orbits x* [16]. We choose the approach comprehensively 
treated by Brunovsky [17], which is particularly amenable to analysis in the extended time- 
delayed feedback case and for which one can then prove the converse of Lemma 3.1: 

extended time-delayed feedback is stabilising for the same gams for which Brunovksy's 
approach is stabilising the classical feedback control. 


Near-impulse feedback and its parametrised monodromy matrix We pick state feedback control in the 
form of a single large but short impulse. That is, we consider a short time 5 G (0, T) and define 
the linear feedback control 


u 6 (t;y) = A 6 (t)K,[y, where 


As(t) 


1/6 ift| mo d[o,T) ^ [0, 6] 
0 if tlmod[0,T) [0/6]. 


(14) 


where t| mo d[ 0 ,T) is the number t G [0, T) such that (t — t)/T is an integer, and Ko G ]R n is a vector 
of constant control gains. Let us first look at classical feedback u(t) = As (t)K ( { [x* (t) — x(t)] 
(where we assume that we know the periodic orbit x*). Using feedback law (14) the feedback 
controlled system reads 

x(t) = f (x(t). As (t)Kg [x* (t) - x(t)]). (15) 


We define the nonlinear time-T map X(x; 6 , Ko) as the solution at time T (the period of the 
periodic orbit x*) of (15) when starting from x at time 0 (including the dependence on param¬ 
eters 6 and Ko as additional arguments of X). Then, for small deviations yo from x*(0), the 
map X(-; 6 , Ko) has the form X(x*(0) +yo; 6 , Ko) = y (T) + O(||yo|| 2 ), where y satisfies the linear 
differential equation (recall that A(t) = 3 x f(x*(t),0), b(t) = 3 u f(x*(t),0)) 

y(t) = [A(t)-b(t)As(t)K ( j] y(t), y(0)=y 0 , (16) 

and the term O(||yo|| 2 ) is uniformly small (including its derivatives) for all 6 . Let us introduce a 
complex parameter y. into (16), which will become useful later in our consideration of extended 
time-delayed feedback: define for a general complex y. with y| G5 C (with an arbitrary fixed 
C > 0) the linear ODE 

y(t) = [A(t) - pb(t)A 6 (t)Ko] y(t), y(0)=y 0 . (17) 

Denote the monodromy matrix of (17) from t = 0 to t = T by P(y.; 6 , Ko) to keep track of 
its dependence on the parameters 6 G (0,T) and Kq G lR n . Thus, P(y; 6 , Kq ) refers to the 
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same monodromy matrix as P(p), defined by (11), for the special case K(t) = (t)Ko. Then 
P(p; 6 , Kq) satisfies 


P(p;6,Ko)=P 0 e X p(-b(0)Kjp) + 0(5), (18) 

where the error term 0(5) is uniform for |p| ^ C and bounded ||Ko||, including its derivatives 
with respect to all arguments. Hence, we can extend the definition of P(p; 6 , Ko) to 5 = 0: 

P(p;0,K 0 ) = lim P(p; 6 ,K 0 ) = P 0 exp(-b(0)Ko p), 

6—>0 


= Po [I — cr(p)b(0)K ( {] where 

(19) 

fexp(^Kjb( 0))—1 
cr(q) = l K 0 T b(0) 0 

( 20 ) 

[h if Kg b( 0 ) = 0 . 



The limit is uniform for all p with modulus less than C. For li = 0, P is the monodromy matrix 
Po of the uncontrolled system, and, thus, independent of 5 and Ko. 

Approximate spectrum assignment for finitely many Floquet multipliers The control (14) is a simplifi¬ 
cation of the general case of finitely many (at most n) short impulses treated in [17]. Feedback 
of type (14) permits us to assign arbitrary spectrum approximately under the assumption that 
the pair (Po,b(0)) is controllable (recall that, according to the definition of Po in (12), Po is the 
monodromy matrix of the uncontrolled system x(t) = f(x(t), 0 ) along the periodic orbit x*). 
This is a stronger assumption than the assumption made in [17], but it is still a genericity 
assumption. 

Lemma 3.2 (Approximate spectrum assignment for classical state feedback control, simplified from [17]) 

Letr > 0 be arbitrary. If the pair (Pq, b(0)) is controllable (that is, the n x n controlla¬ 
bility matrix [b(0), Pob(0),..., Pj l_ 1 b(0)] is regular), then there exist a 5 max > 0 and a 
vector of control gains Ko £ R n in (14) such that all Floquet multipliers of x* for the 
differential equation (15) have modulus less than r for all 5 £ (0,5 max ), where A 5 is as 
defined in (14). 

Note that the vector Ko can be chosen independent of the 5 £ (0,5 max ), but it may depend on 
the radius r into which one wants to assign the spectrum. This result follows from classical 
linear feedback control theory ([17] proves a more general result). In short, linear feedback 
control theory [17] makes the following argument (thus, proving Lemma 3.2): the linearisation 
of X with respect to its initial condition can be expanded in 6 as 

a x X(x*(0);5,Ko) = P(l;6,K 0 ) = Poexp(-b(0)Kj) + 0(6) 

(where P(-; 5, Ko) was the generalised monodromy matrix defined for (17)). Since det Po is posi¬ 
tive we can for every matrix R with positive determinant find a vector Ko such that spec R = 
spec(Poexp(—b(0)K([)) (using the assumption of controllability; see auxiliary Lemma A.l, 
which is a special case from the more general treatment in [17], and [39] for a Matlab imple¬ 
mentation). Hence, if we choose the spectrum of R inside a circle B r / 2 (0) of radius r/2 around 
0, then the spectrum of 3 x X(x*(0); 6 , Kq ) is also inside B r (0) for sufficiently small 5 > 0. 
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Approximate spectrum for extended time-delayed feedback We fix the control gains Ko such that 
lims-^o P(1; 5, Ko) = Poexp(—b(0)Kg ) has all eigenvalues inside B r (0) for some r G (0,1). 
Consider now again the extended time-delayed feedback control (8)-(10) with the particular 
choice of short impulse linear feedback law (14): 

x(t) = f(x(t),A 6 (t)Ko [x(t)-x(t)]) (21) 

x(t) = (1 — e)x(t — T) + ex(t — T), (22) 


where e G (0,1). 

Lemma 3.3 (Floquet multipliers of extended time-delayed feedback) 

Assume that the matrix Po exp(—b(0)K,{) has all eigenvalues inside the ball B r (0) with 
r < 1. Then, for all sufficiently small e and 6, the periodic orbit x(t) = x(t) = x*(t) 
of system (21), (22) has a simple Floquet multiplier A = 1 and all its other Floquet 
multipliers are inside the unit circle. 


Outline of proof (details are given in Section A.2) Eigenvalues A of the linearisation of (21)-(22) are 
roots of the function 


H(A;e,5) 


det 


AI-P 



£ 

A-(l-e) 



(23) 


Roots of h with a non-small distance from 1 — e are close to the roots of det(AI — P(l; 6, Ko)), 
which are inside the unit circle by assumption. Roots A of h close to 1 — £ with modulus 
greater than 1 — e/2 have the form A = 1 — e + e/k where |k| is bounded away from 0 and 
infinity. The roots k of h(l — £ + e/k; 5, e) are small perturbations of the roots k^q of det (I — 
Po — Po1>(0)Kq ct(k — 1)), where a is as defined in (20). These roots k^q have the form 


K £,0 = 1 + 


27ti£ 

kJmo) 


(24) 


(if Kq b(0) / 0, otherwise, only a single root kq/j = 1 exists). The roots k^q have all modulus 
greather than unity (except for t = 0, which corresponds to the trivial eigenvalue A = 1) such 
that the corresponding roots At of h have modulus smaller than unity. (This ends the proof of 
Lemma 3.3.) □ 


Remark — two types of Floquet multipliers The proof of Lemma 3.3 shows that there are two 
distinct types of roots: those approximating the spectrum assigned by the choice of control 
gains Ko, and those close to 1 — e (called At above). The roots At lie close to the circle of radius 
e/2 around the center 1 — e/2 in the complex plane and have the form 


At 



Kjb(0) — 27tif 
Kj b(0) +27tif 


(£ G Z). 


For i = 0, the expression is exact (giving the simple root at unity), for the others the approxi¬ 
mation is sufficiently accurate for small 5 and £ to ensure that they stay inside the unit circle. 
The illustration in Figure 2 shows the two distinct groups for the Hopf normal form example 
discussed in Section 5. 


10 








Importance of scalar input and trivial Floquet multiplier The proof of Lemma 3.3 hinges on one 
argument that depends on the presence of a trivial Floquet multiplier: we need to find the 
roots Sj of s i —y det(I — Po — Pob(0)K ( |s) and then find solutions k of cr( k — 1 ) = sj for all these 
roots Sj. Since b(0)Kj has rank one we know that s ^ det(I — Po — Pob(0)K ( { s) is a first-order 
polynomial (see Section A.2 for details). The presence of a trivial Floquet multiplier then 
ensures that this first-order polynomial has the root 0. Hence, 0 is its only root, restricting the 
possible location for the k^q to the list in (24). This simple argument would not apply for cases 
where the uncontrolled periodic orbit x* has no trivial Floquet multiplier, or for control with 
non-scalar inputs u, or for control with more than one kick per period. 


4. Autonomous feedback control 


The feedback control constructed in Lemma 3.3 introduces an explicit time dependence into 
the system. The controlled system has the form 

*(t) = f(x(t),A 6 (t)Ko [x(t)-x(t)]) 

x(t) = (1 — e)x(t — T) + ex(t — T), 

where A 5 is time-periodic with period T, but the system still has a Floquet multiplier A = 1. 
The neutrally stable direction corresponding to this Floquet multiplier is a phase shift: if 
(x(t),x(t)) = (x*(t),x*(t)) is a periodic orbit of (25) then so is (x(t),x(t)) = (x*(t + s),x*(t + s)) 
for any s £ R. Hence, the controlled system with the gains K(t) = A§(t)Ko is susceptible 
to arbitrarily small time-dependendent perturbations (say, experimental disturbances): the 
phase s of the stabilised solution may drift until A§ is non-zero at a time s where the gains 
Ko are no longer stabilising. This problem does not occur if, instead of applying the feedback 
Kq [x(t) — x(t)] at a fixed time per period, we apply it in a strip in R n close to a Poincare 
section at x*(0) (as illustrated in Figure 3), putting a factor depending on x(t) in front of 
K,{[x(t) — x(t)]. Specifically, we let the function A 5 not depend explicitly on time t but on a 
function t : IR n 1 —y R, where the argument of t is x(t). Then the common notion of asymptotic 
stability of periodic orbits in autonomous dynamical systems applies. One would then always 



x*(t) 



Figure 3: Illustration of choice for strip and Poincare section where gains should be non-zero. 


apply control near x* (0) despite phase drift. A possible explicit expression for u is 
u(t) = A P/ 6 (x(t))Ko [x(t) — x(t)], where A P/ 5 (x) = J p (x)A 5 (t(x)), 

^1 if |x — x*( 0 )| < p, Tr , x*( 0 ) T 

0 if |x — x* ( 0 )| > 2 p, 


JpW = 


t(x) = 


x*(0) T x*(0) 


[x —x*( 0 )]. 


(26) 

(27) 


p > 0 is a small radius and J p is smooth. In (27), t(x* (t)) = t + 0(t 2 ) for |t| <C 1, and J p restricts 
control to the neighborhood of radius p around x* (0). With u as defined in (26), the right-hand 
side of the now autonomous system 

x(t) = f(x(t),A 6 /P (x(t))Kj[x(t)-x(t)]), x(t) = (1 - e)x(t-T) + ex(t-T) (28) 
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has a right-hand side that depends discontinuously on x(t) (because A§ is discontinuous in its 
argument. Since the general mathematical theory for DDEs coupled to difference equations is 
not well developed, one may replace the discontinous A 5 in (26) with a smooth approximation 
of A§. This does not affect the final result, which we can state as a lemma (see Appendices A.3 
and A.4 for the details of the choice for p and the smoothing of A§ p ): 

Lemma 4.1 (Autonomous stabilisability of periodic orbits with extended time-delayed feedback) 

Assume that the matrix Po exp(—bfO )K ( {), as used in Lemma 3.3, has all eigenvalues 
inside the ball B T (0) with r < 1. Then, for all sufficiently small p, there exist e max > 0 
and 6 max > 0 such that the periodic orbit x(t) = x(t) = x*(t) of system (28) is 
asymptotically exponentially stable for all £ £ ( 0 , £ max ) and 6 £ ( 0 ,5 max ). 


Remark: other arguments for A§ /P In (26) we can replace the argument x(t) of A § /P with x(t), 
x(t — T) or x(t — T) without changing the linearisation in x(t) = x(t) = x*(t). Thus, (28) 
successfully stabilises the periodic orbit x* also with these modifications. 

Robustness We assumed perfect knowledge of the periodic orbit x* and the right-hand side f 
in the construction of Ko and A§ /P . However, we know that stable periodic orbits persist under 
small perturbations. Thus, for gains near Ko and functions close to A§ /P the periodic orbit 
of the controlled system persists. Due to the non-invasive nature of extended time-delayed 
feedback, the periodic orbit of the system with perturbed Ko and A 5 /P is still identical to x*. 

5. Illustrative example: Hopf normal form 

The construction of gains as described in Section 4 has been implemented as a Matlab function 
(publically available at [39], depending on DDE-Biftool [40, 41, 42]). The supplementary 
material demonstrates how one can find stabilising gains for two examples: 

1 . a family of period-two unstable oscillations around the hanging-down position of the 
parametrically excited pendulum, and 

2. the unstable periodic orbits in the subcritical Hopf normal form. 

We discuss example 2 in more detail in this section, because for this example we can prove that 
stabilisation with ETDF is not possible with constant gains and small £. The subcritical Hopf 
bifurcation has also been used commonly in the literature as a benchmark example. Here we 
choose the Hopf normal form with constant speed of rotation (such that in polar coordinates 
the angle 0 satisfies 0 = 1 and all periodic orbits have period 27t). Note that the control 
constructed by Fiedler et al [31] depended on changing rotation and was stabilising only in a 
small neighborhood of the bifurcation. Flunkert & Scholl [32] analysed time-delayed feedback 
control (with £ = 1) of the subcritical Hopf bifurcation completely, but also excluded the case 
of constant rotation and restricted themselves to a small neighbourhood of the bifurcation. 
Thus, even though example 2 is seemingly simple, it shows that the method proposed in the 
paper is able to stabilise periodic orbits that are beyond the approaches previously suggested 
in the literature. Without loss of generality we choose a linear control input b = [1,1] T such 
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that the system with control has the form: 


Xi = pxi — X 2 + Xl [x\ + X 2 ] + u, 
X 2 = Xl + px 2 + x 2 [xf + x|] + u. 


(29) 


where p < 0. This system has for u = 0 an unstable periodic orbit of the form x* (t) = 
[r sint, —r cos t] T with radius r = sj —p and period T = 27t. The monodromy matrix Po for the 
uncontrolled system along the periodic orbit x* equals 


Po 


1 0 
0 exp(—47tp) 


Since the derivative of the right-hand side with respect to the control input equals b(t) = b = 
[1,1] T , the periodic orbit is controllable in time T. In fact, the pair (Po, b) is controllable as 
required for the applicability of Lemma 3.3. Extended time-delayed feedback control, applied 
to a two-dimensional system has the form 


u(t) = Ki (t) [xi (t) - xi (t)] + K 2 (t) [x 2 (t) - x 2 (t)] 

Xj(t) = (1 - e)xj(t-T) + exj(t-T) (j = 1,2). 

We can state two simple corollaries from our general considerations. First, it is impossible to 
stabilise the periodic orbit x* with extended time-delayed feedback using time-independent 
gains Ki and K? for small e: 

Lemma 5.1 (Lack of stabilisability for constant control gains) 

Let p < 0 and let Ki(t) and K?(t) be arbitrary constants (also calling them Ki and 
K?). Then there exists an e max > 0 such that the periodic orbit x* is unstable with the 
extended time-delayed feedback control (30) for all e E (0, £ max ). 


Proof Amann & Hooton [38] proved a general topological restriction on the gains K(t) for 
extended time-delayed feedback control: let K(t) E lR n be arbitrary (continuous), 0 <G [0,1] be 
arbitrary, and let u be of the form 

u = 0K(t) T [x(t) —x(t)]. 

The scalar 0 provides a homotopy from the uncontrolled system (0 = 0) to the controlled 
system (0 = 1). Assume that the trivial Floquet multiplier Ai = 1 of x* is isolated for u = 0 
(which is the case for example (29) with p < 0). Then the Floquet multiplier Ai depends 
smoothly on 0 at least for small 0 and will be real: Ai(0) E IR for 0 < 0 < 1. A necessary 
condition for extended time-delayed feedback with gains K(t) to be stabilising for x* and 
arbitrary e E (0,1) is that A] (0) ^ 0 if the number of Floquet multipliers in {z E C : Re z > 1} is 
odd for 0 = 0. If we denote an adjoint eigenvector for the trivial Floquet multiplier by x* (t) 
(the right eigenvector is x* (t)), this criterion can be simplified to 

Jp x4t) T b(t)K(t) T x*(t)dt < o 
Jo x.*(t) T x*(t) dt 
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where b (t) = 9 u f (x* (t), 0) (this simplifying criterion was formulated in general in [33]). For 
our particular example, we have 


x*(t) = x*(t) = 


r cost 
r sint 


b(t) = 


T = 2tt 


and constant gains K| and K 2 such that the necessary condition of [33, 38] is 


Kr + K 2 ^ 0. (31) 

On the other hand, if Ki + K 2 ^ 0 the Jacobian of (29) with classical linear feedback control 

u = Ki(xi,*(t) —Xi) + K 2 (x 2 / *(t) — x 2 ) = Ki(rsint —xi) + K?(—rcost —x 2 ) (32) 

along x = x* (t) has the trace 

tr 0 x f(x(t), K(x*(t) - x(t)))lx{t)—x*(t) = 2p + 4r 2 - Ki - K 2 = -2p - K a - K 2 . 

Since p < 0, this trace is positive if Ki + K? ^ 0 for all t e [0,2re] such that the classical 
linear feedback control (32) cannot be stabilising for the periodic orbit x* = [r sin t, —r cos t] T . 
Thus, Lemma 3.1 implies that extended time-delayed feedback cannot be stabilising either, for 
sufficiently small e > 0 . 

(This ends the proof of Lemma 5.1.) □ 
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Figure 4: Left: Amplitude r = yj— p and unstable Floquet exponent (equals — 2p) along family of 
periodic orbits. Right: gains Ki (note that Ki < 0 always) and K 2 along family of periodic orbits in Hopf 
normal form (29) with feedback law (33). 

Construction of gains For the periodic control gains AsJtJKg (or the autonomous nonlinear 
gains As /P (x(t))Ko) the gains Ko are constructed such the matrix Poexp(—b( 0 )Kg ) has all 
eigenvalues inside the unit circle (for our illustration we choose the target location at ±i/ 2 ). 
Figure 4 shows the amplitude and unstable Floquet exponent of the periodic orbits and the 
gains obtained in this manner (called Ki and K 2 in Figure 4). Since the pair (Po,b(0)) is not 
controllable at the Hopf point, the gains diverge to infinity for p —» 0. In particular, Po = I for 
p = 0 such that it cannot be linearly controllable with a single input. 
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$ assigned (asymptotic) 
o delay induced (asymptotic) 
« numerically computed 


Figure 5: Numerically computed versus asymptotic spectrum for p = —0.25 (e = 0.04, 6 = T/500, 
p = 0.3, Kj = —0.258, K 2 = 4.786).Left: unit circle in the complex plane. Right: zoom into the circle 
around 1 — e or radius e/2. Computed with DDE-Biftool [40,41,42], see supplementary material and 
[39] for the code. 


Illustration of asymptotics Figure 5 shows how the true Floquet multipliers approximate their 
asymptotic values when using the autonomous time-delayed feedback control (44) with gains 
depending on x: 


u(t) = A 6/P (x(t)] [Ki[xi(t] -xi(t)] + K 2 [x 2 (t) -x 2 (t)]] 

kj(t) = (1 - e)xj(t-T) + exj(t-T) (j = 1,2). ( " * 

For the construction of A§ p we used the construction of t proposed in (26)-(27): 

tW = -p-[x-xo] for x E B p (x 0 ), (34) 

yoho 

where yo = [r,0] T , xo = [0,—r] T and r = \/—p. At the particular parameter value p = —0.25 
shown in Figure 5 the uncontrolled periodic orbit is already strongly unstable: the unstable 
Floquet multiplier equals exp(—47tp) ps 23.141. The gains Ki and K 2 , designed to assign the 
Floquet mulitpliers ±i/2, are also large after division by 6. Hence, the range of 5 and e, for 
which stabilisation is successful is small. The values for 6 and e used in the illustration are 
chosen such that deviations from the asymptotic limit are visible but small. 

6. Conclusion and Outlook 

The paper proves conclusively that there are no restrictions inherent in extended time-delayed 
(state) feedback control if one accepts time-periodic gains, while for constant gains general 
positive results are unlikely (as they are absent for classical feedback control of linear time- 
periodic systems). The particulars of the gain construction presented here, following the 
approach of Brunovksy, are merely for the purpose of proving their existence analytically. 
While the result raises the possibility that more general assignment is feasible (since there 
is a lot of freedom in the choice of general time-periodic K(t)) the techniques for proving 
this may have to be different from those in the paper. The central argument of the paper 
rests on the rank-one nature of the control input, making it easy to locate all roots k of the 
transcendental function k > det[I — Pq exp(b(0)K ( | (k — 1))] (the matrix b(0)K ( | has rank one 
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in our case). The argument also exploits the presence of the trivial Floquet multiplier, thus, 
making the result (even if it is based entirely on linear theory) mostly relevant to the analysis 
of nonlinear systems. 


A. Auxiliary Lemmas and detailed arguments of proofs 

Lemma A.1 

Let (A, b) with A E R nXTl and b E ]R n be controllable, let detA > 0, and let 
(Ai,... A n ) E C n be the spectrum of a real matrix with positive determinant. Then one 
can find a vector K e R n such that spec(Aexp(bK T )) = (Ai,..., A n ). 

The proof is given in [17]. A Matlab implementation of the explicit construction is SpecExpAssign. m 
in [39]. 


A.1. Floquet multipliers for extended time-delayed feedback 


Lemma A.2 (Characteristic equation for Floquet multipliers) 

Let A(t) E lR nxn , b(t), K(t) E R n be T-periodic and let £ be positive. Then the Floquet 
multipliers different from 1 — e of the linear system 


x(t) = A(t)x(t) + b(t)K(t) T [x(t) — x(t)] 
x(t) = (1 — e)x(t — T) + ex(t — T) 


are roots of the function 


h : A i—>• det 


AI — P 




(35) 

(36) 


For p E C the matrix P(p) was defined as the solution x at time T of x = [A(t) — pb(t)K(t) T ]x 
with x(0) = I (that is, P(p) is the monodromy matrix of x = [A(t) — pb(t)K(t) T ]x). 


Proof Let A E C be an eigenvalue of the time-T map M of (35)-(36). Let xo, xo (both [— T, 0] >-)• 
R n ) be the components of an eigenvector corresponding to A, and let x\, xi (also both [—T, 0] i-t 
R n ) be the corresponding components of M[x,x]. Then xi(t) = Axo(t) and, by definition of 
the time-T map M, xi(t) = (1 — e)xo(t) + exo(t). Flence, Axo(t) = (1 — e)xo(t) + exo(t), which 
implies (since A / 1 — e) 

£ 

Xo(t) = -- T -rX 0 (t) 

A — (1 — e) 

for t E [—T, 0]. Using this relation, we can solve (35) on the interval [—'T, 0] as 


*o(t) 


A(t)xo(t) +b(t)K(t) T 



x o(t) 


with boundary condition xo(0) = Axo(—T). By definition of the monodromy matrix P this is 
equivalent to 


'■( 1 -rqbij) x «(-T) - Axo(-T). 
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which has a non-trivial solution xo(—T) if and only if the function h in Lemma A.2 is non-zero. 
(This ends the proof of Lemma A.2.) □ 

A.2. Proof of Lemma 3.3 

The characteristic function h(A; e, 6 ) = det[AI — P(1 — e/(A — (1 — e)); 5, Ko)], defined in (23), has 
a root A = 1 for all small 6 and all e G (0,1): a nullvector of I — P(0; 6 , Ko) is x* (0) (corresponding 
to a linearised phase shift). 

For A with modulus larger than 1 — e/2 the term e/(A — (1 — e)) has modulus less or equal 
than 2. Let us pick Si > 0 and a Ci G (0,1) (both small) such that the polynomial 

A i-t det(AI — P(1 — p; 5, K 0 )) 

has all roots inside the ball B( r+ 1 )/ 2 (0) C Bi(0) for all 6 G [0,5i] and p with |p| ^ Ci. This 
is possible since A i-p det(AI — P(l;0,Ko)) has all roots inside the ball B r (0) by assumption 
of the Lemma and the limit of P(1 — p; 5, Kq ) for 5 —> 0 was uniform for bounded p (recall 
lim 5 ^.oP(l — p;5,Ko) = Poexp((p — l)b(0)Kg )). Thus, h(A;e, 6 ) cannot have roots A on or 
outside the unit circle for which 




A — (1 — e) 

holds. Hence, for all 5 G [0, all roots of h(-; e, 6 ) on or outside of the unit circle must satisfy 

e 


Ci ^ 


A — (1 — e) 

We introduce the new variable kgC defined via 


C 2. 


(37) 


A = 1 — e + e/k. 


(38) 


Restriction (37) for A is equivalent to the restriction Ci ^ |k| ^ 2 for k. Hence, for all 6 G [0, Si] 
and e G (0,1), every root A of h(-; 5, e) on or outside of the unit circle corresponds to a root k of 


g ( k; S,e) 


h(l — e + e/k; 8, e) 
£ ' 


det 


1-E+- ) I —P(1 


k; 6 , K 0 ) 


with Ci ^ |k| ^ 2. This one-to-one correspondence of roots of h and g is given via relation 
(38) and includes multiplicity of the roots. Relation (38) also implies that |k| ^ 1, because, 
otherwise, |A| < 1. The function g has a limit 


g(K;S,£) —>• g(K;0,0) for (5, £)->•() 

uniformly for k with Ci ^ |k| ^ 2. Hence, the set of roots k of g(-; 5, e) with Ci ^ |k| ^ 2 is a 
small perturbation of the set of roots of g(-; 0,0): 

g(K;0,0) = det [l-P 0 exp (b(0)K ( j [k — 1])] 

= det [I — Pq — Pob(0)Kg ct(k — 1)] , where 
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a(K —1) 


exp (Kq b(0)(K — 1)) — 1 


Kq b(0) 

K — 1 


if Kg b(0) /O 
if Kjb(O) =0. 


The generalised eigenvalue problem for the matrix pair (I — Po, Pob(0)Kg ) with characteristic 
polynomial cr i-r det(I — Po — Pob(0)Kg cr) is regular because det(I — Po) =0 but det(I — 
Po — Pob(0)Kg cr(—1)) is regular (since Poexp(— b(0)Kg ) has all eigenvalues inside the unit 
circle, I — Po exp(— b(0)Kg ) is regular). As Po b(0)K ( { has rank 1 the characteristic polynomial 
corresponding to cr i-r det(I — Po — Pob(0)Kg a) has degree 1. Moreover, its only root equals 0 
(which must be simple due to the regularity of det(I — Po — Pob(0)K ( { a)). Hence, we know that 
g(K;0,0) = 0 if and only if ct(k — 1) =0. 


Case Kg b(0) =0 If Kg b(0) = 0, this implies that the only root k of g(-;0,0) with Ci ^ |k| ^ 2 
equals unity. Hence, also for sufficiently small e and 5, the only root k with Ci ^ |k| ^ 2 of 
g(-; 5; e) equals unity (since g(l; 5; e) = 0). 


CaseKgb(O) f 0 If Kgb(0) / 0, wehavethat g(i<;0,0) = 0 if and only if exp(Kjb(0)(K — 1)) = 
1 such that the roots are 

27tif 

K,, ° = 1 + Kjb(0) (f£Z) 

The first part of the subscript, £, numbers the roots, the second part of the subscript, 0, indicates 
that 5 = e = 0. Thus, the roots k^q of g( ;0,0) have a modulus 


|k£,oI 



47t 2 £ 2 

(Kg T b(0)) 2 


such that only the roots with index 

\/3 

fmax ^ ^ fmax/ where £ max = — |Kgb(0)|, 

Z7T 

are in the admissible range with Ci ^ |Ke^ol 2. (Hence, both cases, Kjb(0) = 0 and 
Kjb(0) / 0 can be treated equally.) The admissible roots k ( . / 0 of g(-;0,0) are all simple. 
Hence, for sufficiently small 5 and £, g(-; 5, e) will have roots K£ for |£| ^ £ max that are small 
perturbations of k^q, and these roots Kg are the only roots of g(-; 6, e) with modulus in [Ci,2]. 
Since we know that g(l; 5, e) =0, we know that ko = 1 (hence, for £ = 0 the perturbation is 
zero). Furthermore, for non-zero £ with £ C £ m ax/ the modulus of k^o is greater than 1. Hence, 
the perturbed roots k? also have modulus greater then 1 for sufficiently small £ and 6, and 
non-zero |£| ^ £ max . 

Consequently, by relation (38), the only roots of h(-; 6, £) that could be on or outside the unit 
circle are 

A e = 1 - £ + — where |£| ^ £ max . 

Kf 

However, these roots are simple and satisfy Ao = 1 and |Af| < 1 for non-zero £, since | K£ | > 1 
for non-zero £. (This ends the proof of Lemma 3.3.) □ 
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A.3. Details of construction for autonomous feedback gains — Regularisation of the short 
impulse A§(t) 


To avoid discontinous dependence of the right-hand side on the solution, we first regularise 
the discontinuity of the time-dependent gain K(t). Define for 5 G (0, y/T + 1/16 —1/4) (such 
that 26 2 + 6 < T) the regularised version of A§: 


As(t) = < 


1/6 

0 


6 m 


^ + ^ 2 — t| mod [ 0 ; 
62 


h" llmod[0,T) ^ [0/ 

if t|mod[0,T) ^ [6 + 6 2 , T — 6 2 ], 

if t|mod[0,T) ^ (6, 6 + 6 2 ) 


1 ( flmod[0,T) — T + 6 2 \ 

l« m l-*- ) 


if tlmod[0,T) ^ (f ^ /T), 


(39) 


where m : R i —> [0,1] is an arbitrary smooth monotone increasing function with m(s) = 0 for 
s < 0 and m(s) = 1 for s ^ 1. When using A§ as defined in (39) instead of (6) to define the 
linear (now approximately) short-impulse feedback law 

u 6 (t;y) = A 6 (t)Ko"y (40) 


the nonlinear time-T map is still linearisable. Denoting the monodromy matrix of the linear 
system (recall that A(t) = 3 x f(x*(t),0), b(t) = 3 u f(x*(t),0) and using definition (39) for As) 

y(t) = [A(t) - |xb(t)A 6 (t)Ko] y(t), y(0) = y 0 

again by P(y; 6, Ko) then the monodromy matrix of the smoothed system still satisfies (identical 
to (18)) 

P(h; 6, K 0 ) = P 0 exp(—b(0)KQ y.) + 0(6), (41) 

where the error term 0(6) is uniform for bounded y G C and Ko G R n . Hence, we can replace 
the discontinuous definition (6) for As(t) by (39) in (21), and Lemma 3.3 still applies to the 
modified (regularised) system. 


A.4. State-dependent gains K(x(t)) 


Consider a sufficiently small radius p > 0 such that the equation x* (0) T [x — x* (t)] = 0 has a 
unique solution t G R close to 0 for all x G f> 2 P (x* (0)) C R n , thus defining implicitly a smooth 
function t p : 


t p : R n i-> 1R, t p (x) 


root t of x*(0) T [x — x*(t)] =0 if x G B 2p (x*(0)), 

arbitrary such that t p is smooth otherwise. 


(42) 


The function t, defined in (27), is approximately equal to t p along the periodic orbit x* and near 
x*(0): t p (x*(t)) —t(x*(t)) = t —t(x*(t)) = 0(t 2 ). We consider also a regularised indicator 
function for the neighbourhood of x* (0) 


fl if x G B p (x*(0)), 

J p :R n i-^R, J p = < 0 if x ^ B 2p (x*(0)), 

[arbitrary such that J p is smooth ifxG B 2 p(x*(0)) \ B p (x*(0)). 
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and combine t p and J p with A§ as defined in (39) to the smooth globally defined function 

A 6 ,p : IR n i-t IR, A 6/P (x) = J p (x)A 6 (t p (x)) 

When applying A 5 /P to x* (t) the result is identical to the timed impulse A 5 for small 6 : if 
||x* (t) — x* ( 0 ) || < p for all t e [— 6 2 ,5 + 5 2 ] then 

A 6 ,p(x*(t)) = A§ (t) 

for all t e IR. Consequently, the system with extended time-delayed feedback and state- 
dependent gains 

x(t) “ f(x(t),A 6 /P (x(t))Kj[x(t)-x(t)]) (43) 

x(t) = (1 — e)x(t — T) + ex(t — T), (44) 

which is now autonomous with a smooth right-hand side, has for sufficiently small p and 
5 + 5 2 < p exactly the same linearisation along the periodic orbit x(t) = x(t) = x*(t) as 
system (21), (22). The derivative of A§ /P with respect to its argument is multiplied by 0 if 
x(t) = x(t) = x*(t) for all times t in the term As /P (x(t))KQ[x(t) — x(t)] in (43). 

The time reconstruction function t p as defined in (42) satisfies t p (x* (t)) = t as long as x* (t) e 
B 2 p(x*( 0 )). The definition of t p as proposed in (27) in the main text is a 0(t 2 ) perturbation 
of (42) for t of order 6 . Thus, continuity of the Flqouet multipliers implies that the perturbed 
version of t p preserves the stability of the periodic orbit x*. 
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